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PACS 82 . 37 . - j - Single molecules kinetics 
PACS 82 . 35 . Lr - Physical properties of polymers 

PACS 05. 40. -a - Fluctuation phenomena, random processes, and Brownian motion 

Abstract -The effect of self-generated tension in the backbone of a bottle-brush (BB) macro- 
molecule, adsorbed on an attractive surface, is studied by means of Molecular Dynamics simu- 
lations of a coarse-grained bead-spring model in the good solvent regime. The BB-molecule is 
modeled as a backbone chain of L beads, connected by breakable bonds and with side chains, 
tethered pairwise to each monomer of the backbone. Our investigation is focused on several key 
questions that determine the bond scission mechanism and the ensuing degradation kinetics: how 
are frequency of bond scission and self-induced tension distributed along the BB-backbone at 
different grafting density erg of the side chains? How does tension / depend on the length of the 
side chains A'^, and on the strength of surface adhesion Es? We examine the monomer density 
distribution profiles across the BB-backbone at different and relate it to adsorption-induced 
morphological changes of the macromolecule whereby side chains partially desorb while the re- 
maining chains spread better on the surface. Our simulation data are found to be in qualitative 
agreement with experimental results and recent theoretical predictions. Yet we demonstrate that 
the interval of parameter values where these predictions hold is limited in A^. Thus, at high val- 
ues of Es, too long side chains mutually block each other and freeze effectively the bottle- brush 
molecule. 



Introduction. — One of the most outstanding chal- polymers with side chains of poly(n— butyl acrylate) of 
^tenges in modern material sciences is the design and syn- length N ~ 140, self-induced upon adsorption on a solid 
k> thesis of "smart" macromolecules with stress-activated surface (mica), was reported pl [T5HT7] . An experimental 
5_] functions [Tl[5] . During the last decade one observes thus method for control and manipulation of the bond-cleavage 
5^ a rapidly escalating interest in the field of novel poly- in bottle-brush backbones was also proposed [T7]. Thus, 
mer mechanochemistry which, in contrast to the tradi- a selectivity of bond breakage can be achieved by tuning 
tional (non-selective) one, allows to control bond tension the molecular size of such macromolecules which makes it 
on molecular length scales [TJ|3H5]. Examples related to possible to fabricate the brush so as to focus tension in the 
these advances enable, for instance, rupture of specific middle of the molecule. The increase of the bond tension 
chemical bonds [SMZ] , steering the course of chemical reac- in these macromolecules is induced by the steric repulsion 
tivity (BHTTj , changing color of materials [H] or mapping of the side chains as they tend to maximize the number of 
their stress distribution at molecular level |13lfT4] . contacts with the substrate in order to gain energy. This 
In a series of recent experiments, a strong enhancement tension, which depends on grafting density Og, on the side 
of the tension in the (typically, polymethacrylate with de- chain length N, and on the strength of substrate attraction 
gree of polymerization L ~ 3600) backbone of bottle-brush Cs , effectively lowers the energy barrier for bond scission. 
As observed in experiments, self-induced build up of ten- 
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sion proves sufficient to instantly sever covalent bond in 
the backbone. 

The possibihty for breaking strong covalent bonds is also 
an interesting problem from the standpoint of fundamen- 
tal physics. Amplification of tension in branched polymers 
has been considered theoretically by Panyukov and collab- 
orators in several recent works [TM^ by means of scaling 
theory and Self Consistent Field techniques. Numerous 
possible regimes of brush-molecule behavior in terms of 
N, (jg and eg have been outlined [19]. It was argued that 
polymers with branched morphology, physically adsorbed 
on an attractive plane, allow focusing of the side-chain 
tension on the backbone whereby at given temperature 
T the tension in the backbone becomes proportional to 
the length of the side chain, / w fsN. Here /s denotes 
the maximum tension in the side chains, /s « ksT/b, 
with fee being the Boltzmann constant, and b - the Kuhn 
length (or. the monomer diameter for absolutely fiexible 
chains). However, a comprehensive understanding of co- 
valent bond breaking in adsorbed branched polymers still 
has to be reached. Many of the detailed theoretical pre- 
dictions can hardly be measured directly experimentally. 
It remains unclear how exactly tension builds up with 
growing adhesion to the solid surface and impact exercise 
side chain length TV and grafting density Ug on this phe- 
nomenon. Much insight in this respect can be provided by 
computer experiments, yet there have been only few such 
studies [T51[T^ which didn't focus on these questions. 

In an earlier molecular dynamics (MD) simulation study 
of the scission kinetics it was found that mean lifetime of 
a bond becomes more than an order of magnitude shorter 
when the bottle-brush molecule spreads on an adhesive 
surface [18]. These simulation results [18] indicated also 
that the probability distribution for bond scission along 
the backbone of a bottle-brush is sensitive to the grafting 
density (and, therefore, to the degree of steric repulsion) 
of the side chains. The shape of this distribution resembles 
the experimentally established one |15| only for weaker re- 
pulsion when the side chains are short and do not mutually 
block one another. 

In this letter we report our MD results on the first di- 
rect measurement of the adsorption-induced tension / in 
the backbone of a bottle-brush macromolccule whereby we 
examine the influence of varying parameters: grafting den- 
sity (Jg, adsorption strength es/ksT as well as the length 
of side chains N on the resulting increase of /. 

The model. — We consider a three-dimensional 
coarse-grained model of a bottle-brush macromolccule 
which consists L monomers in the backbone connected by 
bonds. Moreover, two side chains of length TV are grafted 
to every a~^-th repeatable unit of the backbone (except 
for the terminal beads of the backbone where there are 
three side chains anchored) . In this way a grafting density 
CTg, which gives the number of side chains pairs per unit 
length, is defined. Thus, the total number of monomers 
in M the brush molecule is M ^ L + 2N[{L - l)ag + 2]. 




Fig. 1 Fig. 2 

Fig. 1: Starting configuration of a bottle-brush molecule (a 
"centipede") with L = 13 (backbone) and A'^ = 3 (side chain), 
so that for a grafting density ag — 1, the total number of 
segments M = 97, and for ag — 1/4, one has M = 43. 

Fig. 2: A snapshot of a thermalized "centipede" with L — 20 
backbone monomers (blue) and 42 side chains (red) of length 
N = 4. Side chains which are too strongly squeezed by the 
neighbors (especially when the backbone bends) are seen to 
get off the substrate in order to minimize free energy. 

The bonded interaction of each segment in the brush with 
its neighbors is described by the frequently used Kremer- 
Grest potential, y«<^(r) = V''™'^ + V^'^^, with the so- 
called 'finitely-extensible nonlinear elastic' (FENE) poten- 
tial: 

V-^-^r) = -O.bkrl In [1 - (r/ro)^], (1) 

and the Weeks-Chandler- Anderson (WCA) (i.e., the 
shifted and truncated repulsive branch of the Lennard- 
Jones potential) given by: 

^wcA(^) = 4e [(a/r)i2 - (a/r)^ + 1/4] e{2^^^a - r) (2) 

In Eq. (HJ fc = 30, tq = 1.5 whereas in Eq. ([2]) one has 
9{x) = or 1 for a; < or X > 0, and e = 1, cr = 
1. The potential V'^'^{r) has a minimum at bond length 
''bond ~ 0.96. Thus, the bonded interaction, V^'^{r), makes 
the bonds of the brush in our model unbreakable which 
is convenient when measuring bond tension is our main 
concern. In the course of the study, however, we have 
also analyzed the case when backbone bonds may undergo 
thermal scission. The bonded potential for the backbone 
monomers is then replaced by V^{r) = V" + V^'^^, with 
a Morse potential 

F"(r) = Z?{l-cxp[-a(r-6)]}2 (3) 

which was recently used in simulations exploring break- 
age of polymer chains [^M^ and bottle-brushes [TH]. In 
Eq. ([S]) a is a constant that determines bond elasticity, 
we use here a = 1. The dissociation energy D oi a given 
bond is measured in units of fc^T, where fcs is the Boltz- 
mann constant and T denotes the temperature. We use 
here D = 1. The monomer diameter in the case of " break- 
able" potential is 6 = 2^^^^a w 1.12ct. All non-bonded 
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interactions between monomers are taken into account be 
means of the WCA potential, Eq. ©. In consequence 
the interactions in our model correspond to good solvent 
conditions. 

The substrate in the present investigation is considered 
simply as a structureless adsorbing plane, with a Lennard- 
Jones potential acting with strength eg in the perpendic- 
ular z-direction, V^- {z) = 4es[((T/2:)^^ - (tj/z)^]. In our 
simulations wc consider as a rule the case of strong ad- 
sorption, es/ksT = 4.0 8.0. With these interaction be- 
tween the monomers, and at different degree of adhesion 
to the surface, we observe well equilibrated bottle-brush 
molecules, as shown in Fig. |3l 




Fig. 3: Snapshots of equilibrated bottle-brushes with L = 121 
backbone monomers and 244 side chains of length A'' = 30 
(cTg = 1) displayed for different adsorption strengths: (a) 
es/kBT=0, (b) es/fcflT = 4.0, and (c) e^/fcisT = 8.0. The 
spreading of the BB-molecule on the surface with growing ad- 
hesion, Es/fes > is clearly seen. 

In reality the thermal energy isA:BT«4.11xlO^^^J/K 
at absolute temperature T sa 300 K, so for a typical Kuhn 
length of cr = 1 nm one obtains for the tension / « 4 pN. 
Estimates show [5T] that this tension is too small to sig- 
nificantly change the lifetime t « 10'^° of covalent bonds 
with energy Ri lOOksT. Adsorbed bottle-brush molecules 
on graphite or water/alcohol mixtures with spreading co- 
efficient 5 w 20 mN/m are capable of generating a force of 
approximately 220 pN [TS] which reduces the lifetime to 
10^ s. In our coarse-grained model we consider bonds with 
dissociation energy of the order of lOfcsT. With Kuhn 
lengths « 1 nm and strength of adsorption e^/fesT « 1-^8 
we find T to be 10^ 10'^ time units (t.u.) whereby 1 t.u. 
is « 10^12 s_ 

The bottle-brush dynamics is obtain by solving a 



Langevin equation for the position q„ = 2/„, z„] of 
each bead in the molecule, 

mq„ = F^„ + Fr*-7qn + RnW, {n,...,N) (4) 

which describes the Brownian motion of a set of bonded 
particles. In Eq. ^ we choose the force with j = M, 
or KG, depending on the character of bonded (break- 
able - Eq. [3] or nonbreakable - Eq. [Ij interaction. The 
infiuence of solvent is split into slowly evolving viscous 
force and rapidly fluctuating stochastic force. The ran- 
dom, Gaussian force R„ is related to friction coefficient 
7 by the fiuctuation-dissipation theorem, 7 = 0.25. The 
integration step is 0.002 t.u. Time is measured in units 
of ma'^/D or -^/mcr^/e (where m denotes the mass of 
the beads, m = 1), depending on whether V^'^{r) or 
V^^ is used. We emphasize at this point that in our 
coarse-grained model no explicit solvent particles are in- 
cluded. Our simulation was performed in the weakly 
damped regime of 7 = 0.25 as we modeled the interface 
between solid substrate and air. However, we carried out 
also simulation in the strongly damped regime for 7 = 10. 
No qualitative changes were discovered except an absolute 
overall increase of the rupture times which is natural for 
a more viscous environment. 

The initially created configurations. Fig. [T] are equili- 
brated by MD for a sufficiently long period of time so that 
the mean squared displacement of the polymer center-of- 
mass moves a distance several (3 5) times larger than 
the polymer size. We then start the simulation with a 
well equilibrated conformation - Fig. [2]- and measure the 
tension / generated by the steric repulsion of side chains 
on the covalent bonds that comprise the macromolecule 
backbone. Tension was measured by sampling the length 
of each bond and using derivative of Eq. (P). In sepa- 
rate runs we allow thermal scission of bonds by using the 
second bonded potential given by Eq. [3l In that case we 
sample the frequency of bond breaking with respect to the 
bond position in the backbone and create a Probability 
Distribution (a rupture probability histogram). 

MD results. — An important aspect of the tension- 
induced scission of bonds in a BB-molecule, adsorbed on a 
solid plane is the probability distribution of rupture (scis- 
sion frequency) along the molecule backbone since it de- 
termines the fragmentation kinetics which can be exper- 
imentally analyzed. In order to interpret experimental 
results [16j[T7] by means of scaling theory [15j[T9], one as- 
sumes that tension of the side chains is focused onto the 
backbone at its ends and transmitted through the middle 
part of the backbone. This would then lead to uniform 
scission probability of each bonds along the backbone. In 
Fig.|4]wc demonstrate that the rupture frequency depends, 
in fact, essentially on the side chain length, N, provided 
grafting density ag is held constant (in an earlier paper we 
found that this frequency does not depend on backbone 
length L for sufficiently long BB-molecules [13). While 
the rupture frequency levels off and stays uniform in the 
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Fig. 4: Rupture probability histograms for a polymer backbone 
with L = 61 presented for different values of grafting density 
(Tg and length of side chains N as indicated. Here es/ksT = 
4.0 and the bonded interaction is governed by Morse potential 
given by Eq. ©. 



middle for Cg = 1 and N = 1, cf. upper left panel in 
Fig. m one observes sharp focusing at both ends of the 
chain for N ~ 5. As argued in Ref. [TS], this effect is due 
to mutual blocking of the side chains which are in the mid- 
dle of bottle-brush molecule. In contrast, side chains that 
are closer to the ends of the backbone are much more mo- 
bile and, therefore, pull on the backbone more frequently. 
Thus, it is not only the tension value but also the rate at 
which this tension is applied that determine eventually the 
distribution of scission rates along the backbone. However, 
for N = 5 and ag = 1/2 the distribution changes qualita- 
tively, and attains a shape as for cTg = 1.0, = 1 again. 
Moreover, for smaller density of grafting, a-g = 1/6, one 
finds "comb-like" distribution which suggests that bonds 
adjacent to the tethered backbone monomers break then 
five times more frequently. 

The distribution of tension / for different values of 
the governing parameters Cs, (Jg and N, is shown in 
Fig. [5^, b, c. Evidently, / reaches a plateau away from 
chain ends. From Fig. [5ji, the mean value (/) is found to 
grow steadily with adsorption strength, eg (for ag = 1), 
and also with grafting density, <7g (for eg/ksT = 8). 
Our simulation results yield a Gaussian distribution of / 
around its mean value (/} (not shown here). The tension 
is seen to focus on the middle of the chain with increas- 
ing side chain length - Fig. [St. One should note again 
the comb-like variation of / over the successive bonds at 
ag = 1/6, i.e. exactly where the side chains are grafted. 

In Fig. |6] we show one of the main results of this study 
- the variation of tension in the bonds along the bottle- 
brush backbone with growing length N of the side chains 
whereby both the adsorption strength, as well as the 
grafting density ag are varied. In all cases / is displayed 
in dimensionless units, taken as the ratio to mean bond 
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Fig. 5: The tension profile along the brush backbones made 
of L = 61 particles displayed as indicated for different: (a) ad- 
sorption strengths es/fcsT, (b) grafting densities (jg, (c) lengths 
of side chains N . (d) Variation of the mean tension (/)/(/o) (in 
dimensionless units) generated in the backbone of bottle-brush 
with L = 121 as a function of adsorption strength talksT and 
grafting density (jg (inset). 



tension of a non-adsorbed BB-molecule /q. Evidently, the 
theoretical predictions |19j for a linear scaling relationship, 
f (X N, appears only for rather short side chains, iV < 10 
at the largest density ag = 1 and adsorption strength 
eg = 8.0 while for weaker adhesion, eg = 4.0, one finds 
/ cx as predicted [12], yet only for 10 < A < 60. 

For loosely grafted side chains, ag < 0.5 and weaker ad- 
sorption, eg = 2.0, the tension becomes progressively in- 
sensitive to A^, saturating at a constant value at shorter 
and shorter length of the side chains. At ag = 1/6 the 
tension / remains unaffected by adsorption whatsoever 
for any N. 

Generally, one observes progressive deviations from a 
power-law relationship for longer side chains N > 70 even 
at the strongest adhesion. We believe that this is due to 
a progressive immobilization of the BB-molecule for long 
side chains which spread nearly parallel to one another, cf. 
Fig. [21 and mutually block each others mobility at large 
values of eg, effectively "freezing" the macromolecule. 
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Fig. 6: (a)-(c) Mean (dimensionless) tension (/) in a bottle- 
brush backbone composed of L = 121 particles is plotted as a 
function of the number of beads in the side chains A'^ for dif- 
ferent adsorption strength talksT (as indicated), and grafting 
density: (a) Og = 1, (b) Og — 1/2, (c) ag = 1/6. Dashed hues 
represent linear fits with different slope. 

Finally, we show the density profile of all BB-mononiers 
in direction perpendicular to the backbone for several 
strengths of the adsorption potential - Fig. [71 The shape 
of such cross-section profiles is important because it has 
been measured experimentally by AFM for poly(n— butyl 
acrylate) side chains on mica (for strong adsorption and 
(Tg = 1), and has also been reproduced by Self- Consistent 
Field theory [TH]. A close inspection of the cross-section 
profiles, shown in Fig. [71 reveals qualitatively the same 
shapes as predicted theoretically [19] : (i) for weak adsorp- 
tion, es/ksT = 2.0, one observes a nearly triangular "tent- 
like" profile; (ii) for moderate adhesion, es/ksT = 4.0, a 
shoulder of monolayer height is seen to form while most 
of the monomers still stay in the cap, and (iii), at the 
strongest attraction, es/ksT — 8.0, one finds a well- 
expressed monolayer and a narrow sharp cap (a cusp). 
Instead of a cusp right above the molecule backbone, how- 
ever, the AFM experiment yields a well-defined parabolic 
maximum [12]. This is most probably due to the defor- 
mation of the BB-cap by the AFM cantilever (about 20% 
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Fig. 7: Monomer density profiles measured perpendicular to 
the backbone of bottle-brushes with L = 121 and A'^ = 30 
displayed for different adsorption strengths as indicated in the 
legend. The inset presents average number of chains above the 
substrate for a brush with L = 61 as a function side chain 
length A'^ for different sets of es/ksT and (jg (triangles stand 
for es/ksT^l and crg = l/2) 

of its height [19]) which thus modifies and distorts its true 
morphology. In fact, the sharp cusp-like cap is formed by 
partitioning of side chains into adsorbed and non-adsorbed 
ones at high grafting density - Fig. [5]d, c. Indeed, some of 
the side chains may relieve the deformation in the strongly 
adsorbed monolayer and minimize the BB free energy by 
lifting off the substrate. In the inset to Fig. [T] one may see 
the variation in the number of such off-chains with length 
N at different degrees of adhesion. For strong adsorp- 
tion, es/ksT — 8.0, the number of off-chains decreases 
with N since the energetic cost of desorption becomes too 
high. One can also see that for ag = 1/2, = 1.0 the 
number of off-chains is close to zero since no crowding of 
monomers on the surface takes place. It should be clear 
that these 'dangling' off-chains bring about the sharp cap 
in the cross-section density profile, observed in Fig. [7] 

Concluding remarks. — In this Letter we report on 
simulation results from our study of bottle-brush macro- 
molecules adsorbed on a structureless attractive surface 
and examine conditions which govern the build-up of bond 
tension and the ensuing frequency of bond scission distri- 
bution along the molecule backbone. Results are com- 
pared to those derived from experiments and theoretical 
considerations. Our findings can be summarized as fol- 
lows: 

• The variation of bond tension / with growing length 
of the side chains A'', adhesion strength eg, and graft- 
ing density Gg comply with some important theoret- 
ical predictions jl9j albeit in a rather narrow inter- 
val. Especially, a scaling law / (x /sA is observed for 
strong adsorption and very short side chains A^ < 10 
only. At moderate adsorption one also finds / cic 
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fsN^^^ in agreement with theory [12], yet for N < 70. 
Further increase of N slows down the tension build- 
up and eventually makes / insensitive to N as the 
side chains mutually block each other and the brush 
molecule gets progressively "frozen". 

• The distribution of tension and the scission frequency 
along the bonds of the bottle-brush backbone change 
strongly with and dg. Different combinations of 
these parameters can induce uniform, comb-like, or 
strongly focused tension distribution in adsorbed BB 
polymers. 

• The cross-section monomer density profile of ad- 
sorbed bottle-brush molecules agrees with experimen- 
tal observations and theoretical predictions of the 
shape at different regimes of adsorption. We find, 
however, the existence of sharp cusp-like peak in the 
profiles which manifests side chain partitioning into 
adsorbed and off-surface chains at strong adhesions 
and large grafting density. 
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